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We investigate the properties of a deterministic walk, whose locomotion rule is always to travel 
to the nearest site. Initially the sites are randomly distributed in a closed rectangular (A/L x L) 
landscape and, once reached, they become unavailable for future visits. As expected, the walker 
step lengths present characteristic scales in one (L — > 0) and two (A/L ~ L) dimensions. However, 
we find scale invariance for an intermediate geometry, when the landscape is a thin strip-like region. 
This result is induced geometrically by a dynamical trapping mechanism, leading to a power law 
distribution for the step lengths. The relevance of our findings in broader contexts - of both 
deterministic and random walks - is also briefly discussed. 

PACS numbers: 89.75.Fb, 89. 75. Da, 05.40.Fb, 05.50.+q, 05.40.-a 



I. INTRODUCTION 

A large number of phenomena in physics, ecology, 
chemistry, economics, etc Q], H, H| , characterized by scale 
invariant distributions, are in many situations associated 
with Levy walks and Levy flights fl, 5, 6|. Furthermore, 
when related to diffusion mechanisms, these types of sys- 
tems present mean square displacements that, for large 
enough times Q, scale as t a with a > 1. 

In such contexts, some of the relevant challenges are to 
determine: (a) if there are global driving forces underly- 
ing the superdiffusive features; (b) how they can emerge; 
and (c) how they are linked to other properties such 
as self-similarity, fractality, noise with f~@ power spec- 
trum and intermittent bursts behavior, ubiquitous in Na- 
ture ||. In fact, many studies address such general ques- 
tions under different perspectives. For instance, one idea 
points to the concept of self-organized criticality [§] . The 
so called spatiotemporal complex systems evolve through 
a series of avalanches towards critical states, which pos- 
sess scale invariance and long range correlations. These 
hierarchical "paths" are unavoidable due to the character 
of the dynamics and are observed in many problems [T3] • 
For deterministic chaotic systems, on the other hand, the 
above queries may be associated either to dynamical frac- 
tional kinetics [2] or to phase space strange non-chaotic 
attractors and attracting sets of particular geometrical 
partitions (see [ll| and refs. therein). 

In the realm of stochastic processes, especially random 
walks, there are different direct causes for superdiffusion 
and power-law tailed decay. To exemplify just a few of 
them, we mention: (i) the evolution governed by frac- 
tional Fokker-Planck equations [l[; (ii) dichotomic sys- 
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terns for which two time scales, microscopic and macro- 
scopic [HJ, can be identified; (iii) the existence of cor- 
relations in the variance of the physically relevant quan- 
tities fl3j : and (iv) random multiplicative processes in 
the presence of a boundary constraint (a repealing bar- 
rier) [3. 

Moreover, it may happen that a Levy or 
power law distribution for the dynamical variables of a 
random process may be a natural way to lead to cer- 
tain outcomes such as: the diversity of species in evo- 
lutionary ecology (l5j ; efficiency optimization in random 
search, e.g., animal for aging 1 1611 in continuous Euclidean 
spaces [13, [H, EH, HH l22l |23H and targets search in 
discrete lattice environments [241 ]; and to avoid the ex- 
tinction edge in scenarios of low availability of energetic 
resources [25| . 

There is a much less studied class of problems known as 
deterministic walks [26|, [2?], [H| . As in the usual stochas- 
tic case, they describe the movement of a walker in a cer- 
tain medium, which can or cannot have a random char- 
acter. However, the rule of locomotion is always taken 
from some purely deterministic model, rather than from a 
probability distribution [28| . Deterministic walks usually 
present the technical difficulties common to nonlinear dy- 
namical s yste ms pfjj and can give rise to superdiffusive 
processes [28( . In fact, they belong to a new class of mod- 
els known as local optimization problems, such as the 
traveling tourist (26[. In contrast to the previously dis- 
cussed examples of purely random walks, it seems that 
for deterministic walks there are no general guidelines 
indicating when the evolution would generate power-law 
distributions for the dynamical variables. 

In the present contribution we study the previous gen- 
eral questions for a specific deterministic walk. We re- 
visit a recently proposed model [29| , in which the walker 
moves in straight lines from site to site, following a "go to 
the closest target site" rule. The sites are randomly dis- 
tributed in a 2D region. As already pointed out in [29| . 
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FIG. 1: Schematics of the search space where the small 
squares represent the randomly distributed target sites. 
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FIG. 2: io as function of L/£ . For L/i > 1 (2D regime), 
£o is the constant 1/y/N, whereas for L/lo <C 1 (ID), io goes 
as 1/(NL). The crossover takes place for L/lo around unity. 



for certain very particular parameter conditions, this 
type of dynamics surprisingly exhibits power law distri- 
bution of step lengths. Here we reveal the mechanisms 
leading to such behavior, not analyzed in (29j . showing 
that the crossover is due to a trapping effect associated 
to particular spatial configurations of the landscape. The 
onset of this phenomenon resembles a critical point in 
thermodynamics, even though there is no real phase tran- 
sition in the system. 

The paper is organized as follows. In Sec. II we pro- 
pose the model. Simulations are presented in Sec. III. 
In Sec. IV we discuss and interpret our findings. Final 
remarks and conclusion arc drawn in Sec. V. 



II. THE MODEL 

We consider a deterministic walk model that was orig- 
inally presented in Ref. [2i| to describe the locomotion 
of spider monkeys during foraging [20]. We define a 
rectangular region of area A and length L\ = L and 
L2 = A/L along the vertical (y-axis) and horizontal (x- 
axis) directions, respectively. Within this domain, a total 
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FIG. 3: The numerically calculated distributions of distances 
between closest neighbor sites for three values of L/£o (see 
main text), corresponding to the: 2D (cross), ID (triangle) 
and crossover region (diamond) cases. The respective analyt- 
ical distributions fitting 2D and ID are plotted as solid lines. 
Note that the intermediate (crossover) is close to the 2D case. 



of N point targets are initially distributed at random. 
The configuration of the search region is schematically 
represented in Fig. 1. In all the simulations we set 
TV = 2.5 x 10 7 and A = 1. The dynamics is given by 
the two following simple rules: 

• Once at a certain target site, the walker moves 
straight to the closest available site. 

• The walker does not come back to any previously 
visited site - the search is destructive, i.e., the total 
number of sites decreases as they are found along 
the walk. 

Let us define the characteristic length 

io = 2rf (1) 

with 

7l = ^J2 d ^ ( 2 ) 

n 

where d n is the distance between the target n and its 
closest neighbor. As L can be taken in the interval [0,1], 
we have two limiting cases. When L = 0(1), the process 
takes place in a 2D space and £0 = y/l/N. On the other 
hand, as L — > the domain is ID and £q = 1/(LN). The 
crossover between these two regimes is found by varying 
L. Fig. 2 displays the numerically calculated Iq as a func- 
tion of L/£q. The two limiting behaviors are clearly seen 
and separated by a crossover emerging around L/£q rj 1. 
In the following, we will use L/£q as the main parameter 
of the model. 

In Fig. 3 we display the distribution D(d/d) of the sep- 
aration distances d n for the three situations: the 2D limit 
with L/£ w 4978.56 (L = 1 and £ = 2.00861 x 10~ 4 ), 
the ID-limit with L/lo ~ 9.99474 x 10~ 3 (L = 2x 10~ 5 
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FIG. 4: Normalized step length distributions for the same pa- 
rameters as in Fig. 3. (a) The triangles (cross) represents the 
ID (2D) limit. The curves show the fits 1.3 exp[-0.92 £/£ ] 
(dashed) and 3.1 (£/£o)~ 5 ' s (continuous), (b) The intermedi- 
ate case (diamond). Here the fit is 0.23 (£/£o)~ 2 ' 2 . 



FIG. 5: (a) The numerical average step length £ (in units 
of £q) taken by the walker during the search as function of 
L/Iq. The region where £ presents a peak, indicated by a 
rectangle, is shown in detail in (b). The continuous curves 
are just guides for the eye. 



and £q = 2.00105 x 10 3 ), and in the crossover region with 
L/£ « 4.21598 (L = 8.82 x 10~ 4 and £ = 2.09204 x 
10 -4 ). The expected distributions, Poisson D(d/d) = 
exp[—d/d] and the standard Wcibull (i.e., a weighted 
Gaussian) D(d/d) — (n/2) d/d exp[— nd 2 / (Ad 2 )], respec- 
tively in the ID and 2D cases, are recovered. 



III. RESULTS 

We simulate the walks by iterating our two rules for 
various values of L/£q. Each simulation runs until the 
walker visits 10 5 targets, out of the initial N = 2.5 x 10 7 . 
The curves are obtained by averaging over 10 3 runs. At 
t = 0, the walker is located on a site in the vicinity of the 
center of the searching environment. 

A first important quantity is the distribution P(£/£q) 
of the reduced distance £/£q traveled by the walker be- 
tween two consecutive targets sites. In Fig. 4(a) we 
show the results corresponding to the ID and 2D cases 
of Fig. 3. For ID, the curve follows the expected ex- 
ponential Poisson distribution. For 2D, P(1/£q) differs 
markedly from the standard Weibull (i.e., a weighted 
Gaussian) distribution of nearest distances D(d/d). The 



curve is broader, but can be well fitted by a rapidly 
decaying inverse power-law with exponent close to 5.3, 
whose general behavior is actually Gaussian, driven by 
the Central Limit Theorem with converging second mo- 
ment Q. However, for the example in the crossover 
region, L/£q sa 4.21598, P clearly exhibits a very long 
tail, as shown in Fig. 4 (b). There is a small but non- 
negligible probability of long walks. In this case, we find 
numerically that P ~ {£/£ )~^ with fx w 2.2 (ji w 2.15 
by considering only the interval 10 < £/£$ < 10 4 ). Thus, 
the distribution has a power-law behavior with a diverg- 
ing second moment, similar to Levy processes. 

A second relevant quantity is the reduced average step 
length £/£o, displayed in Fig. 5 (a) as a function of L/£q. 
As shown in more detail in Fig. 5 (b), we have a peak 
for £/£q in the crossover region, with the maximum cor- 
responding to L/£q w 1.99 « 2. 

From the Fig. 5 we are lead to think that the step 
lengths in the crossover region are indeed larger than 
those for the ID and 2D limits. In fact, we find this is 
true within the interval 2 < L/Iq < 30, where the dis- 
tribution P(£/£q) can fairly be written as (£/^o) _M with 
2 < fj, < 3. For some values of L/Iq, we list in Table I the 
corresponding power-law exponents fj,. Some step lengths 
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FIG. 6: The distribution P(£/£ ) fitted as {l/to)-". The pa- 
rameters are: L/£o = 2.37957 and /i = 2.3155 (open triangle); 
L/£ = 4.21598 and \i = 2.22267 (full circle); L/£ = 9.74390 
and fi = 2.40229 (open circle); L/£ = 22.1615 and fi = 
2.65706 (full square); and L/£ = 38.2288 and p = 2.96385 
(open square). 
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FIG. 7: Angular distribution of the turning angles between 
consecutive steps. The triangles (ID), crosses (2D) and dia- 
monds (crossover region) correspond to the same cases of Fig. 
4. The dotted curve (for the triangles) is just a guide for the 
eye. 



TABLE I: For some values of 2 < L/£ < 30, the fitted /x's 
for P(£/£ ) written as {£/£ y. 
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2.37957 


2.31550 


3.17214 


2.28575 


4.21598 


2.22267 


5.58647 


2.28323 


7.38735 


2.33486 


9.74390 


2.40229 


12.8287 


2.48336 


16.8655 


2.52824 


22.1615 


2.65706 


29.1083 


2.78089 



distributions are shown in Fig. 6. For comparison we also 
plot the example of Fig. 4 (b). 

IV. DISCUSSION 

To understand the above results we turn to the dy- 
namics of the deterministic search process. In the ID 
limit (L — > 0) the random walker tends to follow an al- 
most straight line, with only a few changes in direction, 
mostly occurring during the first steps. On the other 
hand, the 2D limit (L — > 1) is characterized by a much 
larger available space in both directions. Although the 
destruction of previously visited sites makes the walker 
tend to move forwards with higher probability, there is a 
finite fraction of large turning angles along the walk. To 
quantify these features we present in Fig. 7 a normalized 
plot of the angular distribution T(8) of angles between 
two consecutive steps corresponding to the examples of 
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FIG. 8: In the crossover regime, the search space is a narrow 
strip. In such case, a walker moving in average to the right 
may perform many steps in the opposite direction, returning 
near the rightmost visited site through a very long step. 



Fig. 4. In the ID limit, the distribution is very peaked 
at small angle values, indicating that the walker rarely 
deviates from a certain direction (either left or right, de- 
fined shortly after a few initial steps). In the 2D limit 
there is a bias toward the forward direction, nevertheless 
larger turning angles are also likely to happen. 

It is, however, interesting to notice that Fig. 7 alone 
is not sufficient to explain our findings, since the 2D and 
crossover cases present similar T(6). As the search space 
shrinks in one direction, we pass through a crossover re- 
gion from 2D to ID. The singular behavior in this inter- 
mediate regime can be explained in terms of a recurrent 
feature observed in our simulations (see also [1^). In 
this case the walker moves in average towards a given 
direction, say to the right, in a narrow strip-like space. 
However, sometimes it turns to an anti-parallel path to 
visit sites left behind. After some time, the walker ends 
up in a region depleted of targets, which may be far away 
from the rightmost point reached by the trajectory. To 
return to the region rich in unvisited targets located to 
the right, the walker then needs to make a long jump, 
as depicted in Fig. 8. Such mechanism is illustrated in 
Fig. 9, that shows a space-time graph of a simulated 
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FIG. 9: Number of steps vs. the horizontal projection of the 
walker position, in the 2D (a), crossover region (b), and ID 
(c) cases of Fig. 4. Notice the crucial ultra-long steps at the 
crossover geometry (b). 



trajectory in the three different regimes of Fig. 4. 

This dynamical process is particularly sensitive to the 
distance between the two horizontal borders or equiv- 
alently, to the values of L/£q. In fact, for very small 
L/£o (ID) there are no anti-parallel paths, whereas for 
large L/£q (2D) the extra vertical direction often pro- 
vides closer sites than those reached by big jumps across 
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FIG. 10: (a) The numerical drift velocity along x, 
{\x — xo\/n) /£q, as a function of L/£q. The rectangle indi- 
cates part of the crossover region where there is an inflection 
point for the drift, (b) Blow up of the marked region in (a), 
showing a local minimum around L/£o = 1.7094. The contin- 
uous curves are just guides for the eye. 



depleted regions. None of these two aspects, the direc- 
tional bias in the ID case and the extra dimension pro- 
viding many "escape" paths, are present in the crossover 
region. Note that the broad power-law distributions for 
the length of the steps £ (Fig. 4 (b)) are observed when 
the probability of large turning angles is high. In the 
crossover region it is higher than in the 2D regime, as 
shown in Fig. 7. 

The above scenario is confirmed by analyzing two 
quantities related to the dynamics of the determinis- 
tic search process. First, wc calculate the normalized 
drift velocity along the horizontal direction x, defined as 
(\x — xo\/n) /£q, where xq is the starting coordinate and 
x is the coordinate at step n. We show in Fig. 10 (a) the 
drift velocity as a function of L/£q. As expected, it van- 
ishes in the 2D limit. Worthwhile noticing, however, is 
the behavior of the curve in the crossover region, Fig. 10 
(b). In particular, it presents a local minimum around 
L/£q = 1.7094. Furthermore, the first local maximum 
after this minimum is at L/£q = 2.0368 « 2, the same 
position for the maximum of £/£q seen in Fig. 5 (b). 

A second relevant quantity is the fraction of visited 
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FIG. 11: The numerical fraction of visited targets, Xi as a 
function of L/lo. The inset shows a local minimum in the 
crossover region. The continuous curves are just guides for 
the eye. 



targets along the walker trajectory, defined by 



Here, M v - ls , is the average number of visited targets in the 
area searched by the forager and Mo is the total number 
of initial targets in that area. The searched area is de- 
fined by the region [x max - x mia ] x [y max - y min ], where 
the "min" and "max" subscripts stand for the minimum 
and maximum values of the coordinates reached by the 
walker during a full run. Thus, x represents a search 
efficiency of the walker. For the ID regime, all the tar- 
gets are found along the way, so x = 1- On the other 
hand, at 2D regime, x takes a constant small value. The 
crossover from one limit to the other, as a function of 
L/£o, is shown in Fig. 11. Again wc observe a local min- 
imum in the crossover region, as seen in the inset of Fig. 
11. 



V. FINAL REMARKS AND CONCLUSION 

Motivated by the results of a previous study (2j|, we 
have investigated in detail a deterministic walk model 
where the destructive search environment can be changed 
from a 2D to a ID geometry by tuning a single control 
parameter, namely, L/£q. The movement of the walker 
is driven by the "go to closest target" rule. Naively, the 
model should lead to a Poisson process, since the initial 
distribution of target sites (which are destroyed once vis- 
ited) is random. We actually find that in both the ID 
and 2D limiting cases, the step lengths distribution have 
finite variance. However, for some intermediate values of 
L/£q, a non-trivial dynamical process with infinite vari- 
ance takes place, combining a large number of relatively 
small steps with rare long steps. It gives rise to a "Levy- 
like" step length distribution, in which the probability 



of large steps is enhanced by characteristic exponents in 
the range 2 < fjb < 3, e.g., /i « 2.2 for the case in Fig. 
4 (b). Furthermore, for values of the control parameter 
in this crossover region, we observe changes in the search 
efficiency and drift velocity. Such findings are interesting 
since they show that power-law distributions can also re- 
sult from a simple short-range dynamics combined with 
a geometrical constraint. 

Finally, we comment on interesting similarities be- 
tween our model and other related problems. First, a 
random foraging model based on Levy strategies for tar- 
get sites with a recovery (or regeneration or refractory) 
delay time [H, QJ] has been recently studied. Once vis- 
ited, a target site become available for a future visit only 
after a finite number of steps (time delay r). It has, 
as limiting regimes, the destructive (r — > oo) and non- 
destructive (t — * 0) random searches. Differently from 
the present case, the walker may either finish a given 
step with no target found or truncate its flight if a site is 
found along its way. In this foraging problem the most 
efficient destructive (non-destructive) searches requires 
(i ~ 1 (/i f=a 2). So, the parameter r makes the crossover 
from one to other limit. Thus, it seems that the presence 
of boundaries in the deterministic walk and the time de- 
lay in the foraging random search problem play similar 
roles in the sense that they determine the characteristic 
exponent for the distribution of the step lengths of the 
respective walkers, governing then the type of dynamics. 

Second, the foregoing results also suggest an analogy 
between our deterministic walker model and thermody- 
namic systems and phase transitions. The observed mo- 
tion in the 2D limit is relatively isotropic, (of course with 
some bias due to "a back step" depletion), whereas the 
motion in ID breaks completely this isotropy. Similarly, 
the behavior of the velocity is ergodic in the 2-D limit but 
non-ergodic in ID limit. Let us now consider more care- 
fully the crossover. Exactly at the point in between the 
two behaviors, we expect the velocity to be marginally 
ergodic, such that the average behavior of the velocity 
inversion becomes log-periodic rather than periodic (2D 
limit) or nonperiodic (ID limit). Log-periodic velocity 
inversions [30( can represent the border between superdif- 
fusive and diffusive regimes. Moreover, this logarithmic 
behavior implies that mean values of £(t)/£o will scale 
geometrically with time t, such that the mean value of £ 
diverges. In other words, the larger the system size (or 
simulation time) , the larger the mean value of £. Except 
for the memory or correlation effects, this is the same 
kind of behavior we find in Levy walks. The maximum 
superdiffusion for Levy walks occurs when the first mo- 
ment of the mean step size diverges, which corresponds 
to an inverse square distribution of £ . Considered from 
this point of view, the results in Figs. 4 (b) and 5 make 
qualitative sense. 

Thus, although the present study indicates that the 
system is going through a crossover between two different 
limits, from the above discussion we cannot completely 
rule out the possibility of a dynamical phase transition. 
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This issue is presently being investigated and will be re- 
ported in the due course. 
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